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Abstract 

We consider the Einstein-Maxwell equations in space-dimension 
n. We point out that the Lindblad-Rodnianski stability proof ap- 
plies to those equations whatever the space-dimension n > 3. In 
even space-time dimension n + 1 > 6 we use the standard conformal 
method on a Minkowski background to give a simple proof that the 
maximal globally hyperbolic development of initial data sets which 
are sufficiently close to the data for Minkowski space-time and which 
are Schwarzschildian outside of a compact set lead to geodesically 
complete space-times, with a complete Scri, with smooth conformal 
structure, and with the gravitational field approaching the Minkowski 
metric along null directions at least as fast as r^"™ 1 )/ 2 . 



1 Introduction 

There is increasing interest in asymptotically flat solutions of Einstein equa- 
tions in higher dimensions, see e.g. [2,10,14,15]. The pioneering work of 
Christodoulou and Klainerman [4] proving the nonlinear stability of four- 
dimensional Minkowski spacetime uses the Bianchi equations, and therefore 
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does not extend to dimensions larger than four in any obvious way. Now, 
global existence on IR n+1 with n > 4 for small initial data of solutions of 
quasi-linear wave equations of the type of Einstein's equations in wave coor- 
dinates has been proved in [16, 20] 1 , see also [3] for odd n > 5, but the analysis 
there assumes fall-off of initial data near spatial infinity incompatible with 
the Einstein constraints 2 . 

In this note we point out that the Lindblad-Rodnianski stability argu- 
ment [21, 22] in space-dimension n = 3 can be repeated for all n > 3 for 
the Einstein-Maxwell system. Thus, Minkowski space-time is indeed stable 
against electro- vacuum non-linear perturbations in all dimensions n + 1 > 4. 

Next, we point out that non-linear electro- vacuum stability, for initial 
data which are Schwarzschildian outside of a compact set, can be proved by 
the standard conformal method on Minkowski space-time for odd n > 5. 
As usual, the method gives detailed information on the asymptotic behavior 
of the gravitational field, not directly available in the Lindblad-Rodnianski 
method. 

It should be mentioned that in vacuum, and in even space-time dimen- 
sions n + 1 > 4, existence of smooth conformal completions has been proved 
in [1] using the Fefferman-Graham obstruction tensor, for initial data which 
are stationary outside of a compact set. The argument there is simpler than 
the Lindblad-Rodnianski method, but less elementary than the standard con- 
formal method presented here. Moreover, the direct conformal method here 
provides more information about the asymptotics of the fields; however, for 
hyperboloidal initial data our conditions are more restrictive. In any case, it 
is not clear whether the argument of [1] generalizes to Einstein-Maxwell equa- 
tions. (Compare [11, 12] for a completely different approach when n = 3.) 

1 Those works build upon [18,19]; however the structure conditions in [18,19] arc not 
compatible with the Einstein equations. 

2 In [16, 20] compactly supported data are considered. In the theorem for general quasi- 
linear systems given in [3] the initial data are in a Sobolev space which requires fall-off 
at infinity faster than r -™- 3 / 2 . This should be compared with a fall-off of g M „ — rj^ not 
faster than r~ n+2 required by the positive energy theorem. 
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2 Nonlinear stability in higher dimensions 



Consider the Einstein-Maxwell equations, in space-time dimension n + 1, 

D^ v = , 1 ' 

with = ^{T^T* ~ \g^F Xp F\ P ) and = d li A v -d v A (l . It is assumed 
throughout that n > 3. 

We suppose that we are given an n-dimensional Riemannian manifold 
g), together with a symmetric tensor K, and initial data (A = Aidx\ E = 
Eidx' 1 ) for the Maxwell field. (Throughout this work, quantities decorated 
with a bar are pull-backs to the initial data surface y.) For the stability 
results we will assume that y = M ra , but this will not be needed for the local 
existence results. We seek a Lorentzian manifold {^M = Kx y , g) with a 
one-form field A, satisfying (2.1), such that g is the pull-back of g, K is the 
extrinsic curvature tensor of {0} x y, while (A, E) are the pull-backs to y 
of the vector potential A^dx^ and of the electric field T^dx^rf ', where n M is 
the field of unit normals to y . 

We assume the following fall-off behaviors as r = \x\ tends to infinity, for 
some constants a > 0, m: 



j (1 + ^)Sij + O^r' 1 -*) , for n = 3 

9ij - 



Vi,j = l,...,n < 



8ij + 0(r-—- a ) , for n > 4 , 

Ai = 0(r-^- a ) , (2-2) 
K l3 = 0{r-^-«) , 
I ^ = 0{r-^- a ) . 

We will, of course, assume that the constraint equations hold: 

R - K)K{ + K\K] = 2^,' + T, s r : ' , 

V/../ 1 n { 1)< K, y I ), K'j , (2.3) 

ViJ™ = . 

where R is the scalar curvature of g, D is the covariant derivative operator 
associated with g while V is the space-time covariant derivative. 

The following result can be proved by a repetition of the arguments in [21, 
22]. We refer the reader to [23,24] for details and for some information on 
the asymptotic behavior of the fields: 
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Theorem 2.1 Let (y = R n ,g,K,A,E), n > 3, be initial data for the 
Einstein-Maxwell equations (2.1) satisfying (2.2) and (2.3), with ADM mass 
m, set 

(2.4) 

Write g = 5 + + hi with 

h o ixir) 2 ^ forn = 3, 
0l] [0 fom>4, 

for a function x £ C°° equal to 1 for r > 3/4 and to for r < 1/2. Define 3 

E* n „(o)= E (ll(i + 1/2+7+|/| vvXlli 2 + ll(i + ^) 1/2+7+|/| v^||| 2 

0<\I\<N n 

(2.5) 

+ ||(1 + r )i/2+^+|/| vv / j 4||2 2 + ||(i + r ) V2+7+I/I W Ie\ || 2 ) . 

There exist constants e > and 7o(^o); with 7o(^o) as e — > 0, swca 
£na£ 7 /or a// initial data satisfying 

yj E Nn „(0) + m < e Q , (2.6) 

/or a certain 7 > 70, £ae Cauchy problem described above has a global solution 
(g,A) defined on M. n+1 , with (R n+1 , g) — geodesically complete. The solution 
is smooth if the initial data are. 

The threshold value N n in (2.4) arises, essentially, from the n-dimensional 
Klainerman-Sobolev inequalities [17]. 

For initial data which are polyhomogeneous, or conformally smooth, at 
i°, one expects that the solutions will have a polyhomogeneous conformal 
completion at . Unfortunately, the information about the asymptotic be- 
havior of the fields obtained in the course of the proof of Theorem 2.1 does 
not establish this. In the remainder of our work we address this question, 
in odd space-dimension n > 5, for a restricted class of initial data. We will 
use a conformal transformation to both give an alternative, simpler proof of 
global existence for small data, and obtain information on the asymptotic 
behavior. 

3 We take this opportunity to note a misprint in the norm in [23, Eq. (9)]. 



N n = 6 + 2 



n + 2 
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3 Cauchy problem for the vacuum Einstein 
equations in wave coordinates 

We first recall some well known facts. We start with the Einstein equations 
without sources, the Einstein-Maxwell equations are considered in Section 7 
below. 

The vacuum Einstein equations in wave coordinates on W n+1 constitute a 
set of quasi-diagonal quasi-linear wave equations for the components of 
the spacetime metric g, which we write symbolically as 



where the right-hand-side is a quadratic form in the first derivatives of the 
g's with coefficients polynomials in the g's and their contravariant associates. 
We take M n as an initial manifold. The spacetime manifold will be a subset 
V of M n+1 . The geometric initial data on M n are a metric g and a symmetric 
2-tensor K. 

We suppose also given on W 1 the lapse N, shift (3 as well as their time 
derivatives, chosen so that the corresponding spacetime metric and its first 
derivatives satisfy on M. n the harmonicity conditions = 0. If the initial 
data g, K satisfy the constraints, a solution of the Einstein equations in 
wave coordinates on V is such that the harmonicity functions := \3 g x^ 
satisfy an homogeneous linear system of wave equations on (V,g), with also 
d t Fv = 0, hence F M = if (V,g) is globally hyperbolic. 

4 Conformal mapping 
4.1 Definition 

To prove our global existence result we use a mapping : x \— > y from the 
future timelike cone with vertex 0, J+ (0), of a Minkowski spacetime, which 
we denote (M.™ +1 ,r) x ), into the past timelike cone with vertex of another 
Minkowski spacetime, (ffi^ +1 , %)■ This map is defined by, with r\ the diagonal 
quadratic form (—1, 1, 1), 





(3.1) 
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dx a dxP 



x; 



0: 



n+1 

-y 



by x a i-> y a : 



x^ ' 



(4.1) 
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It is easy to check that is a bijection from /+ x (0) onto I yjT) (0) with inverse 
0-1 . y _> x by x a := -^C— . (4.2) 

Moreover is a conformal mapping between Minkowski metrics, it holds that 

r] a pdx a dx^ = n- 2 r] a pdy a dyP (4.3) 

where f2 is a function defined on all R™ +1 , given by 

Q := r) af3 y a yP . (4.4) 

This conformal mapping appears to be better adapted to the context of the 
Einstein equations, which involve constraints, than the Penrose transform 
used in [3] for general quasi-linear wave equations. 

4.2 Transformed equations 

We set 

f»v ■= 9^u - r)txv ■ (4.5) 

We consider / := (f^ v ) as a set of scalar functions on The Einstein 

equations in wave coordinates, for the unknowns / := (f^ u ), are then a 
quasi-diagonal set of quasi-linear wave equations of the form 

The general relation between the wave operator on scalar functions in two 
conformal metrics transforms the left-hand-side of (4.6) into the following 
partial differential operator (compare Remark 4.1 below) 

We introduce the following new set of scalar functions on 

/^rrV/o^- 1 , i.e. /^^n-^of 1 . (4.8) 

With this notation, the system (4.7) reads 
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We note that g alS — i] al3 is a rational function of / = {f^ v ) with numerator a 
linear function of / and with denominator bounded away from zero as long 
as g a /3 is non degenerate. Therefore we have 

(^-^of 1 = Q^h a/3 , (4.10) 

with h al3 a rational function of / and tt 11 ? 1 , with denominator bounded away 
from zero as long as g a p is non degenerate. Now, 

(4.11) 

We use the definitions of Q and of the mapping <fi to compute the right-hand- 
side of (4.9). It holds that: 

dVl 

— = 2r jaf) y 13 =: 2y a , (4.12) 

and 

dv a 

A a »:=/-o ( f ) - 1 = n8Z-2y a y,. (4.13) 



We see that is bounded on any bounded set of IR™ +1 . 

Elementary calculus gives for an arbitrary function / on IR™ +1 the identi- 
ties 

^ Q <P~ l = Af U (4.14) 
dx» 0y a K J 

and 

d2 f „ A -i A p 9 r Aa d(f o (f)- 1 ) 1 
9^ ° =A ^{^ 9r }" ( 4 - 15 ) 

Hence, by straightforward calculation, 

»' ^- BS <Z^ +c .?Uip. ( , 16) 



where the coefficients B and C are bounded on any bounded subset of 
They are given by 

<:=V a ^?> (4-17) 



i.e. 



SJJ = fi 2 5f ^) - 2n(y^) + y^?) + 4y Q y ^ (4.18) 
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and 

= -2fi(y A «£ + y,5 a x + y a Vx J + 8y"y A y M ) . 
If we now set / o = Q k f, we find: 



dy c 

and 

dy a dyP 
with 



dy a dy a 



(4.19) 



(4.20) 



dy a dyP 



+ 2A;fi fc - 1 ( 1//3 |4 + 2/ Q |4) + kn k - 2 D a(3 f , (4.21) 



L> a/3 := 4(fc - l)y Q 2//3 + 2r] Q/3 fi . 
The second term on the right-hand-side of (4.9) is 



(4.22) 



O0- 1 = n^-^F^+n^f) [A^n^l+2kyJ)} 



(4.23) 

Now, A°y Q = — ily M , and recalling that fc = ^y^, the right-hand-side of the 
last equation can be rewritten as 



(4.24) 



This shows that this term extends smoothly, as long as the metric i] + Vt~^~ f 
is non degenerate, to a smooth system on M™ +1 if n > 5 and n is odd. 
The first term at the right-hand-side of (4.9) is 



" I M d^W* + /- 



<9y c 



A similar analysis shows that this again extends smoothly if the metric rj + 
f is non degenerate, if n > 5, and if n is odd. 
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Remark 4.1 The particular case of the conformal covariance of the wave 
equation that we have used, the identity (4-7), results by a straightforward 
computation, when k = t± y- , from the obtained identities, together with 

^WW^ " n " * a W^n + l). (4.25, 
We have proved 

Proposition 4.2 The bijection <f> transforms the Einstein equations in wave 
coordinates on I* x (0) into a quasi- diagonal, quasi-linear system on IZ y (0) 

which extends smoothly to as long as the metric rj + / is non 

degenerate, n is odd and satisfies n > 5. 

n— 1 

Proposition 4.2 immediately shows propagation of the decay rate 
near J? for hyperboloidal initial data which, in wave coordinates, admit a 
smooth compactification as described above with this decay rate. Further, 
hyperboloidal initial data in this class which are sufficiently close to the 
Minkowskian ones lead to solutions which are global to the future of a hy- 
perboloidal surface by the usual stability results. In the remainder of this 
paper we will show that Proposition 4.2 can also be used to construct solu- 
tions which are global both to the future and to the past. 



5 The local Cauchy problem in 
5.1 Initial data 

We consider the Cauchy problem on ffi™ +1 for the Einstein equations in wave 
coordinates with initial data (g, K) on a manifold W 1 embedded as a sub- 
manifold 

M x := {x° = 2A} 

of We suppose these initial data satisfy the Einstein constraints, and 

we choose the initial lapse and shift so that the harmonicity conditions are 
everywhere initially zero. Any globally hyperbolic solution (V x C (?) of 

the Einstein equations in wave coordinates taking these initial data is then 
a solution of the full Einstein equations. Such a solution exists if the initial 
data belong to H 1 ^ x H l ° c , s > f . 
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For our global existence theorem we suppose that the initial data coin- 
cide with a time-symmetric Schwarzschild initial data set outside a ball Br x 
of radius r = R x . Large families of such initial data sets have been con- 
structed in [6,8], arbitrarily close to those for Minkowski space-time. (The 
construction of such data there is done for time-symmetric initial data sets, 
but there is no a priori reason known why all such initial data sets should 
have vanishing K. 4 ) We will use this fact as in previous works [1,5,9], but 
here with the Schwarzschild metric in wave coordinates. 

By the general uniqueness theorem for quasi-linear wave equations, the 
solution (V x ,g) coincides with the Schwarzschild spacetime, in wave coordi- 
nates, in the domain of dependence of the Schwarzschild region, M x \B Rx . 



5.2 The (n + 1)— dimensional Schwarzschild metric in 
wave coordinates 

The Schwarzschild metric g m with mass parameter m, in any dimension n > 3 
is in standard coordinates 

9m = -(l-^) + + r 2 d& (5.1) 



where dD? is the round unit metric on S'™ -1 . Introduce a new coordinate 
system x % = r(r)n\ with n l G S n ; the requirement that x^ = (t,x l ) be wave 
coordinates, DgX^ = 0, is equivalent to the equation 

Setting p = 1/f, one obtains an equation with a Fuchsian singularity at 
p = 0: 

^[p 3 -"(l-2mp- 2 )^] = (n-l)p 1 -V. 

The characteristic exponents are —1 and n — 1 so that, after matching a few 
leading coefficients, the standard theory of such equations provides solutions 



4 A trivial example of initial data set with non-zero K is obtained by moving the initial 
data surface in a time-symmetric space-time, Schwarzschildian outside of a compact set. 
This example raises the interesting question: is it true that vacuum, maximal globally hy- 
perbolic space-times which contain a Cauchy surface which is Schwarzschildian at infinity 
necessarily contain a totally geodesic Cauchy surface? 
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with the behavior 



m 



r = r — 



+ 



(n-2)f n - 3 \ 0(f 5 " 2n ), 



m?_ r 3 In r + 0(r 5 lnr), n = 4 



n > 5 



Somewhat surprisingly, we find logarithms of f in an asymptotic expansion 
of r in dimension n — 4. However, for n > 5 there is a complete expansion 
of r — f in terms of inverse powers of f, without any logarithmic terms: if we 
write g m in the coordinates x\ then 



(9m) (US 

with the functions (/ m ) MI/ of the form 

1 1 x 

(fm)fiv — ~)> x:=(x), 



(5.2) 



(5.3) 



with h IJ/l/ (s,w) analytic functions of their arguments near s = 0. In fact, 
there exist functions h 00 (s), h(s), and h(s), analytic near s = 0, such that 



(- 



1 + 



ioo(r 



„n-2 



/i(r x ) /i(r 1 ) x l xi 



dx l dx 3 . 



As M n is spin, we will necessarily have m > for the initial data that we 
consider. 



5.3 Domain of dependence of the Schwarzschild initial 
data 

The boundary T>+(M x \B Rx ) of the future domain of dependence T>+(M x \B Rx ) 
is threaded by null radial outgoing geodesies of the Schwarzschild metric is- 
sued from B Rx , solutions the differential equation 

dt ( Qm,rr \ 2 /_ . \ 

(5.4) 



dr V g mtt 
such that t(R x ) = 2A, i.e. 



x o = t = 2 \ + I f^H.y d r. (5.5) 

R x ^ 9m,tt ' 
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For the Schwarzschild metric we have 

1 



/ 9m,rr \ 
^ 9m,tt ' 



2m 



< 



Am 

r n—2 



for r > R x provided that 



< 1/2 ; 



(5.6) 



(5.7) 



at fixed R x this can be achieved by requiring m to be sufficiently small; 
alternatively at given m we can increase R x . For n > 5 we deduce from (5.5) 
that on T>+(M x \B Rx ) we have 



b>t-r> a 



with 



b := 2A - R x + 



Am 



:^-3)i?r 3 ' 



2A — R x — 



Am 



[n - 3)i?r 3 



(5.1 



(5.9) 



We can always choose A large enough so that a > 0, requiring e.g. 2 A > R x +1 
for R x > 1 (assuming (5.7)), which implies that V + (M x \B Rx ) is interior to 

/,;,(<))• 



6 The global Cauchy problem 

The local Cauchy problem for the system (4.9) on M y := {y° = — ^-}, with 
initial data in H s+ i x H s , s > | in a ball of radius R y > j^, and i] + f\ My a 
non degenerate Lorentzian metric, has a solution / in a neighborhood of 

V y :=I- y (0)n{y°>-±}, (6.1) 

if the initial data are small enough. The Einstein equations in wave coordi- 
nates on have then a solution on 

V x := ^(Vy) = I^M n x° > A + , (6.2) 

V x is the whole future of the hyperboloid x° = A + VA 2 + r 2 , which is the 
image by of the hyperplane y° = — or- 



12 



The initial data on M y are deduced by the mapping from the values on 
the hyperboloid 

H x := {x° = A + Vr 2 + A 2 } 

of the local solution in and its first derivative, if H x is included in the 
domain V x of its existence. 

The hyperboloid H x is the union of two subsets 

Si := H x n {x° - r > a} and S 2 := H x D {x° - r < a} . (6.3) 

The subset S 2 is included in the Schwarzschild spacetime region. On the 
subset Si it holds that 

A + Vr 2 + A 2 > r + a (6.4) 

A simple computation shows that r is bounded on Si if A > R x and m is 
small, because, using the value of a, one finds that: 

R x (2\-R x )+0(m*) 
T - 2{X-R x -0{m)) ' (6 - 5) 

then x° is also bounded on Si. This subset is therefore included in the 
existence domain V x of the local solution with Cauchy data on M x , for small 
enough Cauchy data. 

We deduce from these results that, for small enough Cauchy data (in- 
cluding small m) on M x the domain V x of the solution contains the future of 
x° = 2A, up to and including H x . 

On <f>(H x fl S 2 ) the initial data f 2 for / is deduced from the Schwarzschild 
metric in wave coordinates, which is static, we have 

kit) = [n-^tfmOtf-W = , (6.6) 

with, using the expression of f m 

I y r \y\ \y\ 

is a smooth function of y (since the origin y = does not belong to the 
domain (f>(H x fl S 2 )) if n > 3. A simple calculation shows that the same is 
true of J^> on 0(if x fl S 2 ) as soon as n > 5. 

On ^(.f/x n Si), the initial data /1, are deduced from the values on the 
uniformly spacelike submanifold H x fl Si of the solution in V x and its deriva- 
tive, by the restriction of to a neighborhood of H x C\Si, where is a smooth 
diffeomorphism. 
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6.1 Conclusions 



By general methods used in [13] one can prove that our global solutions 
are causally geodesically complete. On the other hand, completeness of in- 
extendible spatial geodesies can be established by adapting the arguments 
given in [22, Prop. 16.2], details can be found in [24]. We have therefore 
proved the following theorem: 

Theorem 6.1 Let n be odd and n > 5. Let be given on W 1 gravitational 
data, perturbation of the Minkowski data given by sets of functions f^ and 

■^sffiv Suppose that these data satisfy the Einstein constraints and the initial 
harmonicity conditions and coincide with the Schwarzschild data of mass 
m in the wave coordinates of Section 5.2 outside a ball of finite Euclidean 
radius R x . If these functions are small enough in {H s +i x H s ^ (B(R x )) norm, 
s > n/2, then the data admit a geodesically complete Einsteinian development 
(R n+ \g). 

The smoothness of / as a function of y immediately provides a full asymp- 
totic expansion of g in terms of inverse powers of r = \x\. In fact, we obtain 
an asymptotic behavior of the gravitational field as assumed in [14]. 



7 Einstein Maxwell equations 

Let us show that the conformal method extends easily to the electro- vacuum 
case. 



7.1 The equations 

The Einstein equations with electromagnetic sources (2.1) for the unknowns 
/ := g — 7] take in wave coordinates the form 

where the right-hand-side is a quadratic form in df and dA, with coeffi- 
cients polynomials in g and its contravariant associate. In the Lorenz gauge, 
Va^4 a = 0, the Maxwell equations take also in wave coordinates the form 

a afi d " A p (a)( d f dA ) (72) 

9 dx^- FM{9){ dx~'dx~ h (7 ' 2) 
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where the right-hand-side is bilinear in ^ and in |g, with coefficients poly- 
nomials in g and its contravariant associate. In fact, (7.2) reads in detail 

9 aP -^^ = -{d'tDdpAa - g^ a {d^){d a Ap - d p A a ) ; (7.3) 

this uses both the harmonic coordinates condition and the Lorenz gauge con- 
dition. One thus has a system of equations for (/, A) for which the previous 
analysis applies. 

According to Corvino [7], the existence of a large set of non- vacuum 
electro-vacuum initial data which are Schwarzschildian 5 outside a compact 
set is also valid for the Einstein-Maxwell constraints; compare [9]. 

7.2 Conclusions 

We have obtained the following theorem: 

Theorem 7.1 Let n be odd and n > 5. Let be given on WL n Einstein-Maxwell 
data, perturbation of the Minkowski data given by sets of functions / M „ and 

-^ijfuv, Ax and -J^A\, Suppose that these data satisfy the Einstein- Maxwell 
constraints and the initial harmonicity and Lorenz gauge conditions, and 
coincide with the Schwarzschild data with mass m in the wave coordinates of 
Section 5.2 outside a ball of finite Euclidean radius R x . If these functions are 
small enough in (H s+1 x H s ) (B(R X )) norm, s > | , then the data admit a 
geodesically complete electro-vacuum Einsteinian development (R n+1 ,g,A). 

Similarly to Theorem 6.1, one obtains a full asymptotic expansion both 
of the metric and of the Maxwell field in terms of inverse powers of r near 
Scri. 

Acknowledgements: The authors wish to thank the Albert Einstein 
Institute for hospitality during work on this paper. We are grateful to 
H. Lindblad for bibliographical advice. PTC acknowledges useful discussions 
with J. Corvino. 

5 It could appear that a natural generalisation in this context is to consider initial data 
which coincide with those for the (n + l)-dimensional Reissner-Nordstrom metric outside 
of a compact set. However, if the initial surface topology is W l , then the global electric 
charge is necessarily zero, and one is back in the Schwarzschild case. 
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